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Abstract 

The denning relations (triple relations) of n pairs of parafermion operators (j = 1, . . . , n) 
are known to coincide with a set of defining relations for the Lie algebra so(2n + 1) in terms of 
2n generators. With the common Hermiticity conditions, this means that the "parafermions of 
order p" correspond to a finite-dimensional unitary irreducible representation W(p) of so(2n+l), 
with highest weight (|,|,...,|). Although the dimension and character of W(p) is known 
by classical formulas, there is no explicit basis of W(p) available in which the parafermion 
operators have a natural action. In this paper we construct an orthogonal basis for W(p), and 
we present the explicit actions of the parafermion generators on these basis vectors. We use 
group theoretical techniques, in which the u(n) subalgebra of so(2n + 1) plays a crucial role: a 
set of Gelfand-Zetlin patterns of u(n) will be used to label the basis vectors of W(p), and also 
in the explicit action (matrix elements) certain u(n) Clebsch-Gordan coefficients are essential. 



1 Introduction 

It has been known for a long time that Fermi-Dirac and Bose-Einstein statistics do not yield 
all possible descriptions of particle systems [7, 11]. The classical work of Green [7] is nowadays 
considered as the basic underlying mathematical formulation of the problem of generalized quantum 
statistics. In this paper, the notions of Bose oscillators or bosons and Fermi oscillators or fermions 
are generalized to parabosons and parafermions. Parabosons and parafermions are of interest in 
many applications, in particular in quantum field theory [5,21,23] and generalizations of quantum 
statistics (para-statistics) [7,11-13,16]. 

The "creation and annihilation operators" for parabosons or parafermions satisfy certain triple 
relations, which can be considered as their defining relations. The main object to construct is then 
the generalization of the boson or fermion Fock space, i.e. the paraboson or parafermion Fock space. 
Such a Fock space (with unique vacuum vector) is characterized by a certain parameter p, known as 
the order of parastatistics. One way to construct the Fock space of order p was already discussed by 
Green, and is referred to as the so-called Green ansatz [7]. The mathematical equivalent is related 
to finding a proper basis of an irreducible constituent of a p-fold tensor product [11, 16], in fact 
a superalgebra tensor product for parabosons [16,26]. The computational difficulties arising here 
turn out to be very hard, and did not lead to a complete solution of the problem (constructing a 
proper orthonormal basis of the Fock space, with explicit actions of the paraboson or parafermion 
operators on this basis). 

In a recent paper [17] this problem of giving a complete construction of the Fock space of order p 
was solved for the paraboson case. The solution relies on an important observation by Ganchev and 
Palev [6], who showed that the triple relations for n pairs of paraboson operators are the defining 
relations for the orthosymplectic Lie superalgebra osp(l|2n) [14]. As a consequence, the paraboson 
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Fock space of order p is a certain infinite-dimensional unitary irreducible representation (unirrep) 
of osp(l|2n). The construction of this representation was performed in [17], using group theoretical 
techniques and in particular the branching osp(l|2n) D sp(2n) D u(n). The result is a complete 
description of a proper basis (an orthogonal Gelfand-Zetlin-basis) and the explicit action (matrix 
elements) of the paraboson operators in this basis [17]. 

In the present paper we tackle the same problem for the case of parafermions. It is already 
known for a long time that the triple relations for n pairs of parafermion operators are the defining 
relations for the orthogonal Lie algebra so(2n + 1) [15,27]. The parafermion Fock space of order p 
then corresponds to a unirrep of so(2n + 1), namely to the finite-dimensional unirrep W(p) with 
highest weight (|,|,...,|). Contrary to infinite-dimensional unirreps of the Lie superalgebra 
osp(l|2re), where little was known, a lot is known about finite-dimensional unirreps of Lie algebras 
and in particular of so(2n + 1). For example, it is easy to write down the character of W(p) and 
compute its dimension. On the other hand, a proper basis of W(p), in which the action of the 
parafermion operators can be given explicitly, is not known. The celebrated Gelfand-Zetlin (GZ) 
basis for finite-dimensional irreducible representations of the orthogonal Lie groups or Lie algebras, 
given in [9], is in fact not appropriate here. In the GZ-basis of [9] for so(2n +1), the action is 
given for a set of 2n generators, which are however different from the parafermion operators (and 
also different from the Chevalley generators; in fact, the generators in [9] are not root vectors). 
The relations between the parafermion operators and the generators used in [9] or the Chevalley 
generators are highly nonlinear [25], and hence do not lead to a practical and explicit expression of 
the parafermion matrix elements. 

Previous attempts to solve this problem followed essentially Green's ansatz. In terms of the Lie 
algebra so(2n + 1), this amounts to constructing W(p) as one of the irreducible components in the 
p-fold tensor product W{l)® p of W(l), where W(l) is the so(2n + 1) unirrep with highest weight 
(i, |, . . . , g), i.e. W(l) is the usual fermion Fock space of dimension 2 n (the spinor representation 
of so(2n + 1)). In a number of papers, further attention was given to the actual decomposition 
of W{\)® p into its irreducible components [3,10,24]. In [3,24], the "vacuum subspaces" of the 
irreducible components of W(l)® p were considered, and shown to correspond to unirreps of the 
subalgebra u(n). These approaches did not lead to solution for the parafermion matrix elements, 
however. 

In the current paper we present an explicit solution to this problem. Our solution is based upon 
group theoretical techniques: we use, in particular, the branching so(2n + 1) D u(n), an induced 
representation construction, the known Gelfand-Zetlin basis description for u(n) [8], u(n) Clebsch- 
Gordan coefficients, the method of reduced matrix elements, and the computational techniques 
developed in [17]. Our main result is a specific orthogonal basis for W(p) described by means of 
GZ-patterns, and the explicit action (matrix elements) of the parafermion operators in this basis. 

The structure of the paper is as follows. In section [2] we give the defining relations for the 
parafermion operators, and for the parafermion Fock space W(p). In section [3] we relate the 
parafermion operators to a set of generators for the Lie algebra so(2n + 1), and we give a character 
formula and a dimension formula for W(p). Section [J] is the core of the paper: using an induced 
module construction, we find a proper basis for W(p) and we compute the explicit action of the 
parafermion creation and annihilation operators on the basis vectors. An example is given in 
section [5j and we conclude the paper with a summary and some final remarks in section El 
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2 The parafermion Fock space W(p) 

For a system consisting of n fermions, with creation and annihilation operators F- (i = 1, 2, . . . , n), 
the defining relations are written by means of anticommutators: 

{F-,F+} = 5 lJ , { F r,Fr} = {F l + ,F+} = 0. (2.1) 

The n-fermion Fock space is defined as a Hilbert space with vacuum vector |0), with 

<0|0> = 1, *T|0)=0, (if) f = if (i = l,...,n). (2.2) 

This Hilbert space is finite-dimensional and irreducible under the action of the algebra spanned by 
the elements l,F i FT (i = 1, . . . , n), subject to (|2.ip . A set of orthogonal basis vectors of this 
space is given by 

\e u ...,e n ) = (F+) e - ■ ■ ■ (F+) e "|o>, e n € {o, i}. (2.3) 

This Fock space, denoted by W(l), has dimension 2 n . We shall see that it is a certain unirrep of 
the Lie algebra so(2n + 1), with highest weight (5, • • • , 5) (the spinor representation). 

The purpose of this paper is to study Fock representations of a system of n parafermions. For 
such a system, with creation and annihilation operators (i = 1, . . . , n), the defining relations 
are usually given by 

[//, [/ fc + , /f ]] = 2<Wf , [//> lfk> // + ]] = ^tf - %4 + > [//, [/*. /r]] = 0. 

and their conjugates. They can be written in the following unified form: 

[[ff, fll ffl = \{*- V?hlfl 0%fl (2-4) 

where j, k,l £ {1, 2, . . . , n} and i], e, £ S {+, — } (to be interpreted as +1 and —1 in the algebraic 
expressions e — £ and e — 77). 

The parafermion Fock space W(p) is the Hilbert space with unique vacuum vector |0), defined 
by means of (j, k = 1, 2, . . . , n) 

(0|0) = 1, /r| )=0, Uf) ] = fJ, 

[/7,4 + ]|0)=pM0>, (2.5) 

and by irreducibility under the action of the algebra spanned by the elements /• , fj~(j = l,...,n), 
subject to (|2.4p . The parameter p is known as the order of the parafermion system. For p = 1 the 
parafermion Fock space W(p) coincides with the fermion Fock space W(l). In order to understand 
the structure of W(p) for general (integer) p-values, it will be necessary to make the connection 
between the algebraic relations (|2.4p and the defining relations of the Lie algebra so(2ra + 1). 

3 The Lie algebra 50 (2n + 1) and its unirrep W(p) 

The orthogonal Lie algebra so(2n + 1) consists of matrices of the form: 

a b c \ 

d -a* e , (3.1) 
-e* -c* / 
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where a is any (n x n)-matrix, b and d are antisymmetric (n x re)-matrices, and c and e are (n x 1)- 
matrices (i stands for transpose). Denote by eij the matrix with zeros everywhere except a 1 on 
position (where i,j £ {1, . . . , 2n + 1}). Then the Cartan subalgebra f) of so(2n + 1) is spanned 
by the diagonal elements 

hj = Gjj e n _|_j )ri ,+jf (j = 1, . . . , n), (3-2) 

In terms of the dual basis €j of f)* the root vectors and corresponding roots of so(2n + 1) are given 
by: 

Cjffc — e k+n,j+n £j ~ e fci 3 ' ¥^ k = lj • • • ) n ) 
e j,k+n ~ e k,j+n *~~* e j + e fc; j < k = 1, . . . ,U, 

e j+n,k ~ e k+n,j ~ e j ~ e fc> j < & = 1, . . . , 7t, (3-3) 

ej,2n+l — e-2n+l,j+n *~* ^ji j = 1; • • • > n > 

e ra +j,2n+l — e2n+l,j ^ ~ Cjf, J = 1, ... ,71. 

The positive roots are given by A + = {ej (j = 1, . . . , n); ej — e^, 6j + e k (1 < j < k < n)}. Then it 
is easy to verify that the following multiples of some of the root vectors 

fk = V^C&j^n+l - e 2n +i,n+j), fk = v / 2(e 2 n+i,j - e n+ j i2n +i) (k = 1, . . . ,n) (3.4) 



satisfy the defining triple relations of n parafermions (|2.4p . Furthermore, note that all other root 
vectors of so(2n + 1) can be obtained by considering all possible commutators [ff,f^]. This was 
one of the main results of [15,27]: 

Theorem 1 As a Lie algebra defined by generators and relations, 50 (2n + 1) is generated by the 
2n elements A subject to the parafermion relations 



The parafermion operators fj are part of the positive root vectors, and the f- are part of 
the negative root vectors. The generating vector |0) of the parafermion Fock space W(p) satisfies 
(j,k = l,...,n) 

/"|0)=0, [f7,f+}\0)= P 5 jk \0). (3.5) 
Furthermore, it is easy to verify that 

[fr,f+]=-2h j (j = l,...,n), (3.6) 

from which one deduces hj\0) = — | |0). Hence we have the following: 

Corollary 2 The parafermion Fock space W(p) is the unitary irreducible representation ofso(2n + 
1) with lowest weight (— f , — |, . . . , — |). 

A lot is known about finite-dimensional representations of the Lie algebra so(2n + 1) (or of 
any simple Lie algebra). In particular, dimension and character formulas are available for such 
representations. In these formulas, one uses the Weyl group W and the Weyl tool p. For the case 
of so(2n + 1), W = S n x {—1, +l} n , of order n! • 2 n (S n is the symmetric group). And from (|3,3p 
one can verify that the Weyl tool is 

1 ^ _<A2n + l-2j 2n-l 5 3 1 

aeA + j=l 

Weyl's character formula for an irreducible representation V^(A) with highest weight A reads [29] 

charF (A) = e{w)e< h+p ^/ ^ e{w)e w{p \ (3.7) 
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where e(w) = ±1 is the signature of w, and e stands for the formal exponential. Usually, one 
denotes e ej by Xj\ then the character is a (symmetric) polynomial in the variables xi,x%, . . . ,x n . 

In the current case, W(p) is an so(2n + 1) representation with lowest weight (— §, —5, ■ ■ • , — §)■ 
If p is a positive integer, this is the lowest weight of a finite-dimensional irreducible representation 
with highest weight A = (§, |, . . . , |). Note that, with (ei — £2, ... , e n -i — e n , e ra ) as simple root 
system, the Dynkin labels of this highest weight are [0,0, ... , 0,p\. Let us compute the character 
of this representation: splitting the Weyl group into its symmetric group and {— l,+l} n , the 
numerator in (|3.7p becomes 

£ e(w)e w{A+f,) = £ eM u; j JJ(^+(2n+i-2i)/2 _ ^/2-(2n + i-2 i )/2 ) 
tnsw wes„ y=i 

in terms of a determinant, this can be written as 

Y e(w)e w ^ = det (^+^+1-2^,2 _ P /2-( 2 n+i- 2j )/2 

^ y ' Ki,j<n \ 1 1 

new 

= (si • • • x^-f/ 2 -^"- 1 )/ 2 det(xf +2 ^' - xr 1 )- 



char W(p) = (xi • • ■ x n y^ ae * {X ; 2n _ 3 (3. 



Hence, one obtains 

2 uet^x-. - - x- 
det(x 2n J — xj 

The quotient of determinants appearing here has a nice expression in terms of Schur functions 
s A (x), see [20, p. 84, eq. (2')]. It reads: 

charlU(p) = (xi---x n )~ p / 2 s A (xi,...,x n ). (3.9) 

i{\')<p 

Herein, the Schur functions s\(x) are, as usual [20], labeled by a partition A = (Ai, A2, • • •), where 
all Aj are nonnegative integers and Ai > A2 > • • • • The number of nonzero parts Aj is the length of 
A, denoted by As s\{x) = s\(xi, . . . , x n ) with only n variables Xj, only those A with £(X) < n 
appear (since for ^(A) > n one has s\(xi, . . . , x n ) = 0). In (|3.9p . A' stands for the conjugate partition 
of A (a partition whose Young diagram [20] is the transpose of the Young diagram of A). Otherwise 
said, the sum in (|3.9p is over all partitions A whose Young diagram fits inside the (n x p) rectangle. 
In what follows, the expansion (|3.9p will be very relevant. For a more general context in which 
identities such as (|3.9p appear, see [22] (the expression (|3.9p appears there as [22, Theorem 2.3(1)]). 
The origin of (|3.9p goes back to work of Bracken and Green [2], who do not give (|3.9p explicitly, 
but obtain it implicitly by describing the branching so(2n + 1) D u(n) for W(p). 

Apart from the character formula, there is also Weyl's dimension formula [29] for an irreducible 
representation V(A) with highest weight A: 

dim F A =— . 3.10 

Using (|3.3p and (ei\ej) = 5ij, this yields: 

■ -1 , o • , L(n-1)/2J /p+2n-2i-l\ 

61mW(p)= TT P+1 + 2n " J " fc = TT .fT, . (3.11) 

l<j<k<n J i=0 \2n-4n-l) 

In the last expression with binomial coefficients, the upper limit for i is the largest integer not 
exceeding (n — l)/2. 

The character and dimension of W{p) give already a lot of information, but not yet what we 
really want, namely an explicit orthogonal basis of W(p) and the action of f, on these basis vectors. 
In order to obtain this, we shall follow the induced module construction used in [17]. 
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4 The construction of W(p) 

Prom (|3.3p it is already clear that the 2n elements /• are a set of generators for so(2n + 1). In 
fact, all basis elements of so(2n + 1) (all root vectors plus a basis of the Cartan subalgebra) are 
given by 

ff{l<j<n), [/±,/±] (l<j<k<n), [/+,/-] (l<j,k<n). (4.1) 

Some relevant subalgebras of so(2n + 1) are easy to describe by means of these parafermion gener- 
ators fp: 

Proposition 3 A basis for the subalgebra so(2n) of so (2n + 1) is given by the elements 

[tf,f£](l<j<k<n), [/+,/ fc -] (l<j,k<n). (4.2) 

XTie n 2 elements 

[f+,f k ] (l<j,k<n) (4.3) 
are a frasis /or £/ie so(2n) subalgebra u(n). 

In fact, with [/J 1- ,/^] = 2Ejk, the triple relations (|2.4p imply the relations [E^^E^i] = Sj^En — 
dnEp-j, In other words, the elements [fj,f^] form, up to a factor 2, the standard u(n) or gl(n) 
basis elements. 

So the parafermion generators f^ highlight the subalgebra chain so(2n + 1) D so(2n) D u(n). 
We use here the notation u(n), algebraically the same as the general linear Lie algebra flt(n), but 
with the conditions (fj)^ = fj implying that we are dealing with the "compact form" u(n). 

The subalgebra u(n) can be extended to a parabolic subalgebra [30] V of so(2n + 1): 

V = span{[/+, /"] (1 < j, k < n), fr (1 < j < n), [/r, /"] (1 < j < fc < „)}. (4.4) 

Recall that [/~, /^]|0) = pSjk |0), with \fj, ft] = —2hj. This means that the space spanned by |0) 
is a trivial one-dimensional u(n) module C|0) of weight (— |, . . . , — |). Since /j r |0) = 0, the module 
C|0) can be extended to a one-dimensional V module. Now we can define the Verma module or 
the induced so(2n + 1) module W(p): 

W(p) =Indp (2n+1) C|0). (4.5) 

This is an so(2n + 1) representation with lowest weight (— |, . . . , — §). By the Poincare-Birkhoff- 
Witt theorem [4], it is easy to give a basis for W(p): 

(ft)" 1 • • • / 2 + ]) fcl2 ([/i + , / 3 + ]) fc13 • • • ([/n + -i, /„ + ]) fc - 1 '"|0), (4.6) 

Note that W(p) is infinite-dimensional. It is in general not an irreducible representation of so(2n + 
1). Let M(p) be the maximal nontrivial submodule of W(p). Then the simple module (irreducible 
representation), corresponding to the parafermion Fock space, is 

W{p) = W{p)/M{p). (4.7) 

For W(p), it is easy to compute its character. Using the fact that the weight of |0) is (— |, . . . , — f ), 
the basis (|4.6p . and the Cauchy identity [18] 

nnr _ m i_ JT — r J2s^,..., Xn)= ^ ix) , (4.8) 
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where the sum is over all partitions A and s\(x) is the Schur symmetric function, one finds 



char W(p) 



(xi 



-p/2 



( Xl ---x n )-P/ 2 Y,sx(x). (4.9) 

A 



These Schur functions s\(x) are the characters of finite-dimensional u(n) representations. Hence 
this expansion yields the branching to u(n) of the so(2n + 1) representation W{p). This gives an 
elegant possibility to label the basis vectors of W(p). For each irreducible representation of u(n) 
one can use the corresponding Gelfand-Zetlin basis [1,8]. The union of all these GZ basis is then 
the basis for W(p). Thus the new basis of W(p) consists of vectors of the form (the label p is 
dropped from the notation of the vectors) 



\m) = \m) n 



mi n 

"ll,n-l 



\ 



J 





[m] n \ 




\m) n - 1 J 



(4.10) 



The top line of this pattern, also denoted by the n-tuple [m] n , is any partition A (consisting of 
non increasing nonnegative numbers). The remaining n — 1 lines of the pattern will sometimes 
be denoted by |m) n_1 . All rriij in the above GZ-pattern are nonnegative integers, satisfying the 
betweenness conditions 



> rriij > m i+ i )j+1 (1 < % < j < n - 1). 



Since the weight of |0) is (—§,..., — |), the weight of the above vector is determined by 



h k \m) 



k fc-i 
+ ^2 m jk ~^2 m j,k-i ) \m). 
i=i i=l 



(4.11) 



(4.12) 



Now we use the same technique as in [17]. The triple relations (|2.4p yield 

With the identification [f^, fj~) = 2Eij in the standard u(n) basis, this is equivalent to the action 
Eij • efc = (5jfcej. Hence the set (fx, f^i ■ ■ ■ > fn) ls a standard u(n) tensor of rank (1, 0, . . . , 0). So 

ne can attach a 

-6j. Explicitly: 



one can attach a unique GZ-pattern with top line 10 • • • to every fj~, corresponding to the weight 



10 • • • 000 
10---00 

0---0 





(4.13) 



where the pattern consists of j — 1 zero rows at the bottom, and the first n — j + 1 rows are of the 
form 10 • • • 0. The tensor product rule in u(n) reads 



(H n ) ® (10 • • • 0) = ([mftO © (Hl 2 ) © • • • © i[m]l n ) 



(4.14) 



where ([m.] n ) = (mi n ,m2„, . . . ,m nn ) and a subscript ±fc indicates an increase of the A;th label by 



±1: 



(H± fc ) = {mi n , m kn ± 1, . . . , m nn ). 



(4.15) 
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In the right hand side of (|4,14p . only those components which are still partitions (i.e. consisting of 
nondecreasing integers) survive. 

Now a general matrix element of can be written as follows [17]: 



{m'\ff\m) 



\m\ 



rn 



to 



An— 1 



in 



,n-l 



J .l 

10' 
10' 





to 



\m) 

•00 
•0 



n— 1 



rn 



in 



+k 

An— 1 



X ([to] 



to 



(4.16) 



The first factor in the right hand side is a u(n) Clebsch-Gordan coefficient [17,28], the second 
factor is a reduced matrix element. By the tensor product rule, the first line of |m') has to be of 
the form (|4.15p . i.e. [m'] n 



[in 



+k 



for some A:- value. 



The special u(n) Clebsch-Gordan coefficients (CGCs) appearing here are well known. They can 
be found, e.g. in [28]. They are expressed by means of u(n)-u(n — 1) isoscalar factors and u(n — 1) 
CGC's, which on their turn are written by means of u(n — l)-u(n — 2) isoscalar factors and u(n — 2) 
CGC's, etc. The explicit form of the special u(n) CGCs appearing here is given in Appendix A 
of [17]. 

Just as in [17], the main problem is now to find expressions for the reduced matrix elements, 
i.e. for the functions G k ([m] n ), where 

G k {[m] n ) = G k (m ln ,rn 2n ,...,in nn ) = ([m]" fc ||/ + ||[m] n ), 

for arbitrary n-tuples of non increasing nonnegative integers [m] n = (TOi n ,TO2 n , 
case, one can write the explicit actions: 



,m. 



(4.17) 
In that 



ff\m) 



fj M 



k,m' 



£ 

k,m' 



[TO | 

\<m 
[to 

I TO 



in— 1 



k 

A 71-1 



10- --00 
10---0 



10 • • • 00 
10- --0 





TO 



771. 



/\n— 1 



G k ([m} n ) 



[TO 
TO 



+k 

/\n—l 



[m] n 




[mf_ k \ 


to)™ -1 


j G k {[m] n _ k ) 


\m') n ~ l J 



(4.18) 



(4.19) 



The key calculation to determine the unknown functions G k is to start from the following action: 

n n—1 

[fnJn}H = -2h n \m) = (p- 2(Y j m jn — ^2 m j,n-i))\ m )- (4.20) 

j=l 3=1 

Now one can express the left hand side by means of (|4,18p - (|4.19p . using the explicit form of the 
CGCs and isoscalar factors (which are rather simple in the case j = n). The result is a complicated 
system of coupled recurrence relations for the functions G k - It is quite surprising that this system 
of coupled nonlinear recursion equations has a rather simple solution. Using the relevant boundary 
conditions, we have been able to solve this system of recurrence relations. This task would have 
been hardly impossible without the use of Maple. Our main computational result is: 

Proposition 4 The reduced matrix elements G k appearing in the actions of /• on vectors \m) of 
W{p) are given by: 

G k (mi n ,m 2n , ■■■ , m nn ) = 

(£ n (m kn +n — k) + l) U.j^ k =i( m kn - m jn - k + j) 



n 



Ln/2j 
3*\=1 



m kTl 



m 2 j, n -k + 2j)(m kn - m 2 j, n -k + 2j + l) 



(4.21) 
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for k even; and by 



G k (mi n ,m 2n , ■ ■ ■ ,m nn ) = 

( n \ V2 

(p - m kn + k - l){O n (m kn + n-k) + l) ]l"^fc=i ( m kn - m jn -k + j)\ 

r~~ /2l I V / 

l n^fc+i =1 ( m fcn - m 2 j-i,n - + 2j - l)(m fcn - m 2 j-i, n - k + 2j) J 
fork odd. 

Herein £ and O are the even and odd functions defined by 

Sj = 1 if j is even and otherwise, 

Oj = 1 if j is odd and otherwise; (4.23) 

where obviously Oj = 1 — Sj, but it is still convenient to use both notations. Also, note that 
products such as Il^fc=i m eans "the product over all j-values running from 1 to n, but excluding 
j = k" . The notation [a] (resp. \a}) refers to the floor (resp. ceiling) of a, i.e. the largest integer 
not exceeding a (resp. the smallest integer greater than or equal to a). 

To present all the details of this computational result is unrealistic. Just as in [17], the proof 
consist of verifying that all triple relations (|2.4p hold when acting on any vector \m). Each such 
verification leads to an algebraic identity in the n variables of the partition, mi n , . . . , m nn . In such 
computations, there are some intermediate verifications: e.g. the action [fj~, fk~]\ m ) should leave 
the top row of the GZ-pattern \m) invariant (since [f^,f k ] belongs to u(n)). In fact, it must 
give (up to a factor 2) the known action of the standard u(n) matrix elements Ejj. in the classical 
GZ-basis. 

The explicit expressions for the reduced matrix elements give the action of the generators in 
the basis of W(p), for arbitrary p. The structure of the maximal submodule M{p) and hence of 
the irreducible factor module W(p) is revealed by examining when these matrix elements vanish. 
It follows from (|4.2ip - (|4.22p that the only crucial factor is 

(p — mkn + k-1), (k odd). 

In particular for k = 1 this factor is (p — mi n ), and m\ n is the largest integer in the GZ-pattern. 
Starting from the vacuum vector, with a GZ-pattern consisting of all zeros, one can raise the entries 
in the GZ-pattern by applying the operators . However, when m\ n has reached the value p it can 
no longer be increased. As a consequence, all vectors \m) with m\ n > p belong to the submodule 
M{p). This uncovers the structure of W(p): 

Corollary 5 For p a positive integer, the parafermion Fock space W{p) has as a basis the vectors 
\m) with min < V- ^ n other words, the top line of\m) is a partition A with largest part not exceeding 
p, i.e. i(X') <p. As a consequence, it follows from (14.90 that 

char W(p) = (xi ■ ■ ■ x n )~ p/2 ^ s x {xi, . . . , x n ). 

A, i(X')<p 



This last result coincides with (13. 9p . Note that p must be a positive integer in order to have a 
positive inner product (m\m) for all possible patterns. 

Of course, we did not make these hard computations just to find the character of W(p), which 
was already determined by standard techniques. Our main purpose was to find a proper orthonor- 
mal basis for W(p), in which the action of the parafermion operators f^ 1 can be computed explicitly. 
This has now been obtained by means of (|4.18p - (|4.19p and (|4.21|) - (|4.22|) . Let us summarize this 
result, inserting also the explicit CGCs of [17, Appendix A]. 
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Theorem 6 An orthonormal basis for the parafermion Fock space W(p) is given by the vectors 
\m), see (|4.10p . with m\ n < p. Its dimension is given by (|3.1ip . The action of the Cartan algebra 
elements of so (2n + 1) is: 



k k-l 



hk\m) = { -- + ^m jk - J \m), (1 < k < n). (4.24) 

3=1 3=1 

The action of the parafermion operators f^ is given below, first for j = n (simple case) and then 
for general j. We have: 

fn \ m ) = i>2 ( —^-^77 \ ^ ] Gi{mi n ,m 2n ,---,m nn )\m) +in ; (4.25) 

i= l \ llk^i=l\ L kn - Hn) J 

/n M = V ( ™ k=1 ^' n ~ 1 , lm \ : ) Gi{m ln ,...,m in - 1, . . . ,m nn )|m)_ in ; (4.26) 

Herein, lij = mij — i, Gi is determined by (|4.2ip - (|4.22p . and ±i,n attached as a subscript to \m) 
indicates a replacement mi^ n — * mi^ n ± 1. In general, 

n n - x ' (nr = \(h J -i-i h , J -i)\ 1/2 



f+\m)=^2 ^ ...^2s(i n ,i n -i)S(i n -i,i n - 2 )...S(ij + i, 



i n =li n -.=l i,=l V llW=lV.*-H, 



x 



n 



Y[k^i n - r =l(^k,n-r ~ h n - r +l,n-r+l ~ 1) Y[k^i„- r+ i=l(^k,n-r+l ~ h n -r,n—r) 



\\Ylk^i n - r+1 =l(^k,n-r+l h n -T+l,n-r+l)Y\k^in-r=^ k ' n ~ r k n - r ,n-r 1)/ 

x G in (m ln ,m 2n , . . . , m nn ) \m) + i n>n - + i n _ un -i-.... + i j j; (4.27) 

n n-1 j / rP _1 (Z fc ! — /j ■) \ ^ 2 

= X] • ' ' 5 , (^n^n~i)5'(i n _i,i n _2) • • .5(^+1," "•' 



j-j / rifc^j n _ r =l(^,n-r — k n - r+1 ,n-r+l) Y\k^i^ r+1 =l (h,n-r+l ~ k n - r ,n-r + 1)\ 
r=l \n^i„_ r+ i=l(^,n-''+l ~~ ^„_ r+ i,n-r+l + 1) rifc^j n _ r =l Vk,n-r ~ h n - r ,n-r) J 



X 



x G in (mi n , . . .,rrii n>n - 1, . . . ,m nn ) Im)-^;-^.^-!;...;-^-. (4.28) 

Once again, each symbol ±ik, k attached as a subscript to \m) indicates a replacement rrii k ^ — > 
m iktk ± 1, and 

5 Example: W(p) for so (5) 

The above formulas for general n look rather involved, so it is perhaps useful to give an example. 
Let us consider the case n = 2, i.e. the Lie algebra so (5). For any positive integer p, the parafermion 
Fock space W(p) has dimension given by (13. lip . 

AiraWip) = (i±m±3k±H, ( 5.i) 

6 
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The set of orthonormal basis vectors of W(p) is given by all vectors with integer m-patterns of the 
form 

m 12 ,m 2 2 



V > m vi > mn > m 2 2 > 0. 



(5.2) 



Clearly, the number of different m-patterns satisfying these inequalities coincides with (|5.ip . The 
action of the parafermion generators /{*" and f% is given by: 



ft 



ft 



mi2,m 2 2 
mn 



mi2,m 2 2 
mn 



=V (mn ~ m 2 2 + - mtt) 
- \J (mi2 - mn)(m 2 2 + 1) 
^ =V (mu - mn + - mi 2 ) 
+ \/ (mn ~ m 2 2)(m22 + 1) 



mu + l,m 2 2 
mn + 1 

m 12 ,m 22 + l 
mn + 1 

mi2 + l,m 22 
mn 

mi 2 ,m 22 + l 
mn 



(5.3) 



(5.4) 



The action of /{" and / 2 ~ follows from the above expressions using (m'\f~\m) = (m\f^\m'). In this 
case, it is not too difficult to check these actions by hand: one can do this by comparing the action 
of the left hand side of the triple relation (12. 4D with the action of the right hand side of (12. 4p (for 
some values of j, k, I and £, m e). 



6 Summary and conclusion 

In this paper we have given a complete description of the unitary irreducible representations W(p) 
of 5o(2n + 1) with highest weight (|, |, . . . , |), which in particular are of interest in the theory of 
parafermion statistics because these representations are exactly the Fock spaces of n parafermions, 
with p being the order of the parastatistics. Although many aspects of W(p) are known from general 
Lie algebra representation theory, the problem of giving an appropriate orthogonal basis for W(p) 
and the explicit action of the parafermion operators on such a basis was not solved. In the present 
paper we solve this problem using some group theoretical methods and computational techniques. 
A crucial role in our analysis is played by the u(n) subalgebra of 50 (2n + 1), generated by all 
commutators of the parafermion creation and annihilation operators. Taking a certain parabolic 
subalgebra V containing u(n) and a trivial module of V generated from the vacuum, i.e. the lowest 
weight vector of weight (— |, — |, . . . , — f ), an induced module W(p) of so(2n + 1) is constructed. 
The Fock module W(p) is the quotient of this induced module by its maximal submodule M(p). The 
character of the induced module is readily obtained and by a classical result of Cauchy-Littlewood 
the characters can be rewritten as an infinite sum over partitions of Schur symmetric functions. 
This can be reinterpreted as a decomposition of the so(2n + 1) module into an infinite sum of 
finite-dimensional simple u(n) modules labeled by partitions. For each irreducible representation of 
u(n) one can use the corresponding Gelfand-Zetlin basis. The union of all these GZ basis vectors is 
the basis for the induced module W(p). The main calculation is then the action of the parafermion 
operators on this basis. Just as for the case of paraboson operators [17], the collection of the 
n parafermion creation operators is a u(n) tensor operator of rank (1,0, ... ,0). To calculate its 
matrix elements, they are written as a product of certain u(n) Clebsch-Gordan coefficient and a 
reduced matrix element. As the relevant u(n) CGCs are known, the problem is to find the reduced 
matrix elements. Solving a set of recurrence relations for these, leads to their expressions. These 
explicit expressions give not only the action of the generators in the basis of W(p), they also yield 
the structure of the maximal submodule M(p) and hence of the irreducible factor module W(p). 
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This leads to the main result of the paper: an explicit basis of W(p) (consisting of all possible 
GZ-patterns with integer entries at most p) and the explicit action of the generators in this basis. 
As an illustration, the case n = 2 is given in more detail. 

We have considered here the parafermion Fock spaces of order p for a finite degree of freedom 
(n finite). The real interest lies in such quantum systems (parabosons and parafermions) with an 
infinite degree of freedom (n = +00), see e.g. [11,16]. In a forthcoming paper, we hope to report on 
how the results of [17] and of the current paper can be used to construct representations of order 
p for an infinite set of parabosons and parafermions. Also an investigation of representations of 
the "parastatistics algebra" in which both parabosons and parafermions appear [19] should be of 
interest. 
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